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Rarefield flows about hypersonic vehicles entering the 
upper atmosphere or through nozzles expanding into a near 
vacuum may only be simulated accurately with a direct 
simulation Monte Carlo (DSMC) method. Under this Cooperative 
Agreement, the models employed in the DSMC method were 
enhanced, and simulations in support of existing NASA projects 
and missions were performed. 

DSMC models were developed and validated for simulating 
rotational, vibrational, and chemical relaxation in 
high-temperature flows, including effects of guantized 
anharmonic oscillators and temperature dependent relaxation 
rates. State-of-the-art advancements were made in simulating 
coupled vibration/dissociation/recombination for post-shock 
flows. 

Models were also developed to compute vehicle surface 
temperatures directly in the code rather than requiring 
isothermal estimates. These codes were instrumental in 
simulating aerobraking of NASA's Magellan spacecraft during 
orbital maneuvers to assess heat transfer and aerodynamic 
properties of the delicate satellite. 
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NASA also depended on simulations of entry of the Galileo 
probe into the atmosphere of Jupiter to provide drag and flow 
field information essential for accurate interpretation of an 
on-board experiment. 

Finally, the codes have been used extensively to simulate 
expanding nozzle flows in low— power thrusters in support of 
propulsion activities at NASA Lewis Research Center. Detailed 
comparisons between continuum calculations and DSMC results 
helped to quantify the limitations of continuum CFD codes in 
rarefied applications. 


During the first four years of this research program, the 
work was performed under the direction of Dr. Iain Boyd. Dr. 
Brian Haas joined the team in January, 1991, and became sole 
investigator in January, 1993. The results obtained during 
the phase when Dr. Boyd was Principal Investigator were 
presented in periodic research reports dated 11/2/89, 5/15/90, 
11/28/90, 4/5/91, 3/6/92, and 1/4/93 as summary report. The 

scientific papers published in the open literature resulting 
from this phase of the program are included in the listing 
below. 
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The first periodic report with Dr. Haas as Principal 
Investigator was submitted on 12 October, 1993. His research 
from 1/1/93 through 1/31/94 included a detailed study of the 
rates of rotational and vibrational relaxation in the DSMC 
method to resolve numerous inconsistencies and differing 
interpretations appearing in the literature. This study 
developed a consensus among many researchers in the field and 
was reported in Ref. (15). 

A second research area of Dr. Haas concerned modeling of 
gas-surface interactions. Rather than require prescribed 
temperature estimates for spacecraft surfaces, as is typically 
done in DSMC methods, a new technique was developed which 
couples the dynamic surface heat transfer characteristics into 
the DSMC flow simulation code to compute surface temperatures 
directly. This model, as applied to thin planar bodies such 
as solar panels, was described in Ref. (13) . 

Application of the DSMC method to problems of practical 
interest requires a trade-off between solution accuracy and 
computational expense and limitations. A parametric study was 
performed to assess the accuracy penalties associated with 
simulations of varying grid resolution and flow domain size, 
and was reported in Ref. (16). 
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The DSMC code was applied to two different spacecraft 


entry studies. First, the pitch, yaw, and roll aerodynamics 
of the Magellan spacecraft during entry into the Venus 
atmosphere at off— design attitudes was studied and reported in 
Ref. (14). Secondly, a study of entry of the Galileo probe 
into the atmosphere of Jupiter was conducted to assess vehicle 
drag in support of the on-board Atmosphere Structure 
Experiment during the most rarefied portion of the trajectory. 
This material was submitted in the form of an abstract for the 
AIAA Thermophysics Conference slated for June, 1994 (Ref. (17) 
— see Appendix) . 


The following papers and publications resulted from the 
research activities conducted under Cooperative Agreement 
NCC2-582 : 


(1) I.D. Boyd, " Assessment of Chemical Nonequilibrium in 
Rarefied Hypersonic Flow," AIAA-paper 90-0145. 

(2) I.D. Boyd, P.F. Penko and L.M. Carney, " Efficient Monte 
Carlo Simulation of Rarefied Flow in a Small Nozzle ," 
AIAA-paper 90-1693. 

(3) B.L. Haas and I.D. Boyd, "Vibrationally-Favored 
Dissociation Applicable to a Particle Simulation , " 
AIAA-paper 91-0774. 

(4) B.L. Haas and J.D. McDonald, " Validation of Chemistry 
Models Employed in a Particle Simulation Method , " 
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AIAA-paper 91-1367. 


(5) i.D. Boyd and T. Gokcen, "Evaluation of Thermochemical 
Models for Particle and Continuum Simulations of 
Hypersonic Flow," AIAA 27th Thermophysics Conf., 
Nashville, TN, 1992, AIAA-paper 92-2971. 

(6) I.D. Boyd and E.E. Whiting, "Decoupled Predictions of 
Radiative Heating in Air Using a Particle Simulation 
Method," AIAA 23rd Plasmadynamics & Lasers Conf., 
Nashville, TN, 1992, AIAA-Paper 92-2954. 

(7) I.D. Boyd, "Analysis of Vibration-Dissociation-Recombina- 

tion Processes Behind Strong Shock Waves of Nitrogen," 
Phys. Fluids A 4 (1), 1992, 178. 

(8) B.L. Haas and D. Schmitt, "Simulated Pitch, Yaw, and Roll 
Torques on the Magellan Spacecraft During Aerobraking, " 
Atmospheric Flight Mechanics Conf., 1993. 

(9) B.L. Haas, "Particle Simulation of Satellite Aerobraking 
with Coupled Surface Heat Transfer," 18th Rarefield Gas 
Dynamics Symposium, 1993. 

(10) B.L. Haas, "Flow Resolution and Domain of Influence in 
Rarefield Hypersonic Blunt-Body-Flows , " AIAA 28th Thermo- 
physics Conf., Orlando, FL, 1993; AIAA Paper 93-2806. 

(11) I.D. Boyd and E.E. Whiting, "Decoupled Predictions of 
Radiative Heating in Air Using Particle Simulation 
Method," AIAA Paper 92-2971, Presented at the AIAA 23rd 
Plasmadynamics & Lasers Conference in 6/92. 

(12) B.L. Haas, D. Hash, G . A . Bird, F.E. Lumpkin, III, and 
H . A . Hassan, "Rates of Borgnakke-Larsen Thermal 
Relaxation in Direct Simulation Monte Carlo Methods, " 
Phys. of Fluids A, submitted. 

(13) B.L. Haas, "Models for Dynamic Surface Temperatures 
During Rarified Aeropass Maneuvers, " AIAA 28th Thermo- 
phys. Conference, Orlando, FL, 1993; AIAA-Paper 93-2765; 
presented at the AIAA Thermophysics Conference in 7/93. 

(14) B.L. Haas and D . A . Schmitt, "Simulated Rere field Aerody- 
namics of the Magellan Spacecraft During Aerobraking, " 
AIAA Paper 93-3676; presented at the Atmospheric Flight 
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Mechanics Conference in August, 1993, and submitted to 
J. of Spacecraft and Rockets. 


(15) B.L. Haas, D.B. Hash, G . A . Bird, F.E. Lumpkin, and H.A. 
Hassan, " Rates of Thermal Relaxation in Direct Simulation 
Monte Carlo Methods," Phys. of Fluids A, 1994 (in print) ; 
presented at the Fluid Dynamics meeting of the APS, 

11/93 . 

(16) B.L. Haas M.A. Fallavollita , " Flow Resolution and Domain 

of Influence in Rarefied Hypersonic Blunt-body Flows , " 
AIAA-paper 93-2806; presented at the Thermophysics Conf. 
in 7/93 and submitted to J. Thermophysics and Heat 
Transfer . 
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Simulated Rarefied Entry of the Galileo Probe 
into the Atmosphere of Jupiter 

Brian L. Haas* 

Eloret Institute , Palo Alto, California 94303 
Frank S. Milos f 

NASA Ames Research Center, Moffett Field, California 94035-1000 

Flow properties and aerodynamics are computed with a direct simulation Monte Carlo (DSMC) 
method for rarefied entry of the Galileo Probe into the atmosphere of Jupiter. Accurate predictions of 
vehicle drag coefficients are needed in order to assess atmospheric properties from the onboard Atmo- 
sphere Structure Experiment where highly-sensitive accelerometers will measure the drag force to within 
1 0-6 barr during the initial entry phase at high altitudes. The corresponding flow rarefaction extends from 
the free molecule limit to the near continuum transition regime (Re<1000). Simulation results indicate 
that C D varies from 2. 1 at the free molecule limit down to L6 at Re^ = 1 , 000. Temperatures , densities , 
and internal energies throughout the flow field were also computed at each altitude ranging from 735 km 
to 353 km above the 1 barr level in the Jovian atmosphere. Surface heating and temperatures of the probe 
were computed directly in the DSMC code by assuming radiative equilibrium. Material response was 
re-assessed accurately during entry by accounting for conductivity ; heat capacity, and pyrolysis which 
led to surface material mass efflux several times that of the free stream mass influx. The simulation also 
accounted for the quantum nature of the rotational energy mode of the dominant atmospheric species H 2 
through partial internal excitation in the freestream gas. 


NOMENCLATURE 

A thermal accommodation coefficient 

C D drag coefficient 

D diameter 

gj degeneracy of rotational quantum level j 
j rotational quantum level 

k Boltzmann constant, 1.3805 x 10 -23 J/K 

Kn Knudsen Number 

M Mach number 

q net convective heat flux 

Re Reynolds number based on diameter 
rj ' normalized rotational energy of level j 
T temperature 

t time during entry, starting at 735 km 

Z collision number for internal relaxation 
a VHS exponent of intermolecular potential 

e material radiative emissivity 

Cr rotational degrees of freedom 

9 characteristic mode temperature 

p viscosity (kg/m 2 -s) 

p mass-density (kg/m 3 ) 

cr Stefan-Boltzmann const., 5.67 x 10“ 8 W/(m 2 K 4 ) 

* Research Scientist, Aerothemiodynamics Branch. Member, AIAA. 
Mailing Address: NASA Ames Research Center, M/S 230-2, 
Moffett Field, CA 94035-1000. 

| Research Scientist, Thermal Protection Materials Branch. 

Moffett Field, CA 94035-1000. 

This piper is declared a work of the US. Government and is not 
subject to copyright protection in the United States. 


Subscripts: 

d deep-space value 

ref reference value 

r pertains to the rotational mode 

v pertains to the vibrational mode 

w value at wall or surface 

w value at wall or surface 

oo freestream value 

INTRODUCTION 

Just prior to encountering and orbiting Jupiter, the 
Galileo spacecraft will release a probe which will enter the 
Jovian atmosphere. During the initial aerobraking phase, 
the 45-degree blunted-cone probe will be protected from 
heating by a carbon phenolic shield. Once the probe veloc- 
ity has been reduced from 47.5 km/s to 0.74 km/s, and entry 
heating has diminished, the probe will eject its heat shield 
and deploy a parachute. During descent, the probe will 
make several in situ measurements of atmospheric proper- 
ties and transmit that data to the orbiting spacecraft How- 
ever, an Atmospheric Structure Experiment, 1 similar to that 
employed in the Pioneer Venus mission, will also be on 
the Galileo probe to measure deceleration during the ini- 
tial high-altitude entry phase. This experiment deduces 
atmospheric density, pressure, and temperature from de- 
celeration measurements so long as the vehicle drag co- 
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efficient is known a priori . The instrument is sufficiently 
sensitive to detect any deceleration exceeding 10~ 5 m/s 2 . 
Consequently, meaningful properties can be assessed for 
the Jovian upper atmosphere where the probe encounters 
highly rarefied flow during entry just prior to peak heat- 
ing and ablation of the heat shield. This flow regime is 
bounded by the effective free molecule limit at 750 km al- 
titude (Re^ = 0. 1) and the near-continuum limit at 350 km 
(Re^ = 1, 000). Note that since Jupiter has no identifiable 
surface, altitude is measured relative to the 1.0 barr pressure 
level in the Jovian atmosphere. 

Accuracy of the experiment, however, depends upon the 
accuracy with which the probe drag coefficient is estimated. 
Intrieri 2 conducted a series of experiments in the ballistic 
range facilities at NASA Ames Research Center to mea- 
sure the drag of several blunt-body configurations includ- 
ing the probes from Pioneer Venus and Galileo. Results for 
each were very similar, leading to values near C D = 1 . 1 for 
Re^ > 1,000. However, C D rises significantly (approach- 
ing values near 2.0) for decreasing Re^ below 1,000, al- 
though Intrieri’s experiments were not reliable in the rar- 
efied regime. Furthermore, the experimental surface mate- 
rials, gas species, and flow conditions (Mach number, den- 
sity, etc.) differed significantly from those anticipated for 
the probe entry at a given Reynolds number. 

Due to the lack of sufficient experimental data, the probe 
aerodynamics must be estimated computationally. Unfor- 
tunately, the flow regime is ill-suited to simulations which 
are based upon the continuum Navier-Stokes equations due 
to limitations in the constitutive relations for heat flux and 
shear stress. Instead, highly rarefied flows, for which the ra- 
tio of molecular mean free path to a body dimension is large 
(Knudsen number, Kn>0. 10), are best simulated computa- 
tionally with direct simulation Monte Carlo (DSMC) par- 
ticle methods. 3,4 Here, gas dynamics is modeled directly 
by the motion and interaction of thousands or millions of 
discrete particles. Particles which strike the vehicle may 
reflect back into the flow with velocities and internal ener- 
gies corresponding to full or partial accommodation to the 
surface. The simulation permits accurate assessment of ve- 
hicle aerodynamics and heating along with properties of the 
flow field. 

The objective of the present study was to assess the drag 
coefficient for the Galileo probe during entry from 750 km 
down to 350 km altitude in the Jovian atmosphere. Due to 
uncertainties in the applicable surface thermal accommo- 
dation coefficient A, simulations were repeated using dif- 
ferent values to quantify its effects upon the vehicle entry. 
Surface heating was evaluated to determine the extent of 
pyrolysis of the heat shield through the use of the Charring 
Material Thermal Response and Ablation (CMA) program. 
This code models transient convective heating, radiation, 
in-depth conduction, heat capacity, and the flow of pyroly- 


sis gases through the porous material. Together, these codes 
simulate the entry environment and response of the Galileo 
probe to estimate the appropriate vehicle aerodynamics re- 
quired for this and other aerobraking missions. 

SIMULATION MODELS 

The DSMC code employed in the present study was de- 
veloped by Baganoff and McDonald 4 - 5 and enhanced for 
better application to rarefied aeropass maneuvers. 6 The flow 
field is divided into cubic cartesian cells to facilitate selec- 
tion of colliding particles and sampling of macroscopic flow 
properties. The body geometry is modeled by a compos- 
ite of planar facets in those cells through which the body 
surfaces pass. Each surface facet collects statistics regard- 
ing momentum and energy flux and may assume a surface 
temperature independent of neighboring facets. To simu- 
late entry with the DSMC method, one must first specify 
properties of both gas-gas and gas-surface interactions. 

Gas-Gas Interaction Models 

Molecular interaction is simulated by the Variable Hard- 
Sphere (VHS) model of Bird 7 in which the collision out- 
come corresponds to isotropic scattering, akin to the me- 
chanics of hard sphere interactions. The collision rate, how- 
ever, corresponds to an inverse power-law intermolecular 
potential of exponent a. This parameter must be specified 
between the limits of the Maxwell molecule (a = 4) and 
the Hard sphere (a = co), and may be estimated from the 
known temperature dependence of gas viscosity p as fol- 
lows, 


A* 



( 1 ) 


Transport properties of the Jovian atmosphere, thought to 
be composed by a mixture of 89% H 2 and 11% He, were 
calculated from kinetic theory by Biolsi. 8 That work em- 
ployed sophisticated semiempirical interaction potentials to 
solve the detailed collision integrals for viscosity, thermal 
conductivity, and binary diffusion for the gas mixture. Bi- 
olsi’s results for viscosity are plotted in Fig. 1 and com- 
pared to curves corresponding to the VHS model. The 
VHS parameters which yield best agreement in the fig- 
ure and were employed in the current work are given by 
/iref = 2.24 x 10“ 5 kg/m-s, TW = 1000 K, and a = 10.5. 

Of additional concern in gas-gas interactions is the in- 
elastic exchange of molecular energies due to relaxation of 
the internal energy modes for rotation and vibration. The 
mechanics of these exchange processes are modeled in the 
DSMC code by the methods of Borgnakke and Larsen 9 and 
Haas, et al . 10 These involve partitioning post-collision ther- 
mal energies in a manner which corresponds to equilibrium 
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distributions. The rates of rotational and vibrational re- 
laxation are dictated by collision numbers Z T and Z v , re- 
spectively. Solutions of the master equation using state- te- 
state transition probabilities estimated from quasi-classical 
methods 11 “ 13 have led to rotational relaxation rates in H 2 
described by Z r = 100. Simple application of the Millikan 
and White 14 empirical expression for vibrational relaxation 
rates suggests that Z v is several orders of magnitude larger 
than Z r , leading to exceedingly slow vibrational excitation. 
Furthermore, given the high characteristic temperature of 
vibration (0 V = 6320 K) n for H 2 , it is unlikely that the 
vibrational mode will contribute significantly to the ther- 
mophysics of this flow and was therefore neglected in the 
present work. 

DSMC methods typically model molecular rotation as a 
continuous energy mode with two degrees of freedom, £ f = 
2. However, at freestream conditions, the quantized nature 
of the rotational mode for H 2 leads to an effective num- 
ber of degrees of freedom below two as a result of its high 
characteristic temperature (S r = 85.33 K). 15 As plotted in 
Fig. 2, dependence of ( r upon temperature T is described 
by the rotational partition function and may be expressed as 
follows, 


EJi r i ex P ( — r i) 


Cr=2^ 


£<7;exp(— rj ) 

j 


( 2 ) 


where the normalized energy rj and degeneracy gj per ro- 
tational quantum level j are defined as follows, 16 


rj =j(j + 1)^2 (3) 

9j — 2j + 1. W 

Rather than encorporate a detailed model of the quan- 
tized rotational mode into the code, a simpler adaptation of 
the continuous energy model was employed in the present 
work. First, it was recognized that collisions promoting 
rotational excitation will occur predominantly in the high- 
temperature regions of the flow. Since the rotational de- 
grees of freedom would become fully excited in these re- 
gions, the normal Borgnakke-Larsen model for exchange 
mechanics was employed assuming ( r = 2. However, 
since the freestream temperature was fairly low, incom- 
ing particles were initialized such that the mean rotational 
energy was given by ( r kT /2 where ( r < 2 is evaluated 
from Eq. (2). This prevents the freestream from being 
too enei^etic in rotation, such as would occur if initialized 
with ( r = 2. 


Gas-Surface Interaction Models 

Details regarding the interaction of gas molecules with 
surfaces are not well understood in general, and are depen- 
dent upon several factors including surface roughness, im- 
pact dynamics, molecular potentials, and thermal energies. 
However, simple engineering models often suffice to simu- 
late interaction phenomena. In the present work, a single 
thermal accommodation coefficient A describes the frac- 
tion of particles which accommodate fully to the surface 
versus those which reflect specularly. Full accommoda- 
tion implies diffuse reflection of the particle from the sur- 
face with thermal energy corresponding to the surface wall 
temperature, T w . Rather than assign some temperature to 
the surface, the code couples a simple surface heat transfer 
model into the flow solution to compute T w directly. 17,18 
This model assumes that each surface facet is in radiative 
equilibrium with space at temperature T D = 150 K, leading 
to the energy balance given by 

q — ter (7^ — T* d ) - 0. (5) 

Here, q is the net convected heat flux to each facet account- 
ing for both incident and reflected energy. As will be de- 
scribed later, this model can be enhanced to account for ma- 
terial heat capacity, thermal conductivity through the heat 
shield, and pyrolysis of the surface material. 

Definitive values for A and £ for the interaction of Jo- 
vian atmospheric gases upon carbon phenolic material at 
flight conditions are not readily available. As used in the 
this study, a surface emissivity of £ = 0.85 was suggested by 
Bueche 19 from ground-based experiments, flight data, and 
theoretical predictions, and is at least consistent with esti- 
mates cited elsewhere. 20-23 Appropriate values for A may 
be found with even less certainty than e from limited related 
experimental data 24 and theoretical 25 results. It therefore 
proved necessary to repeat simulations using different val- 
ues of A to assess sensitivity of the results to this parameter. 

SIMULATION RESULTS 

The DSMC code was used to simulated entry of the 
Galileo probe at several points along its trajectory from 
735 km altitude to 353 km. Table 1 lists the simulation 
conditions for each case with atmospheric data taken from 
Orton. Note that time is measured relative to the 735 km tra- 
jectory point. For all cases the velocity was assumed to be 
47,450 m/s and the Knudsen number and Reynolds number 
were based on the probe diameter (1.265 m). The grid res- 
olution employed in the present work is defined in Table 1 
by the size of the probe diameter measured in cell-lengths, 
and was sufficiently fine to yield less than 1% error in drag 
and heating 26 The geometry of the probe is compared to the 
simulation models in Fig. 3. Note that only one quadrant of 
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the probe was simulated, taking advantage of two planes of 
symmetry of the body. 

In general, flows with greater Re^ require greater res- 
olution in order to resolve flow gradients and avoid over- 
predicting drag and heat transfer. However, the required 
size of the flow domain increases with lower Re^ because 
molecules which reflected from the probe surface are capa- 
ble of diffusing far into the flow when collisions are scarce. 
The extent of the upstream diffusion of particles is depicted 
in the plot of flow temperature along the stagnation stream- 
line in Fig. 4. The upstream domain boundary for each sim- 
ulation case was sufficiently far upstream that the transla- 
tional temperature was near its freestream value to prevent 
overprediction of heat transfer and drag. 26 Density profiles 
along the stagnation streamline are plotted in Fig. 5 and also 
exhibit the effects of rarefaction. Note that no clear shock 
structure is observed since the shock is fully merged with 
the body layer. Density rose considerably near the body 
surface due to particle reflection from the relatively cold 
surface. 

For each case, the simuladon employed at least 16 par- 
ticles per cell in the freestream and roughly 3 to 7 mil- 
lion particles total. The code was optimized for vector- 
processing on Cray supercomputers, requiring roughly 
0.6 ^sec/particle/timestep on the Cray C-90 or a total run 
time of 5,000-8,000 CPU seconds depending upon the case. 
Memory requirements ranged from 50 to 200 megawords. 

Results of each case, identified by Re^ , are presented 
in Table 2. The cases were run with .4 = 0.75 and several 
were repeated with A = {0.5, 0.9}. Heating of the nose re- 
gion was highly sensitive to accommodation coefficient A 
while drag was fairly insensitive, particularly for the high- 
est and lowest Re^ cases. Simulated drag coefficients for 
the Galileo probe are plotted in Fig. 6 and compared to the 
experimental ballistic range results of Intrieri for spheres 
and the Pioneer Venus probe. Unfortunately, Intrieri’s re- 
sults for the Galileo probe were all at high angles of attack 
and were of questionable quality at low Re^ . 

In general, drag dropped with increasing Re^ above 10, 
but did not appear to blend well with the high-Re^ exper- 
imental data- However, the experiments were performed at 
lower Mach numbers (roughly M^ = 14.5) in CO 2 com- 
pared with the high Mach numbers expected for Galileo en- 
try in the H 2 -He Jovian atmosphere. Indeed, for comparing 
highly rarefied flows, a suitable parameter for correlating 
the data would be the Knudsen number Kn^ . Re-plotting 
the data in Fig. 7 suggests a smoother transition between 
the simulation results and the experimental work. 

Nonetheless, it was instructive to simulate a few of In- 
trieri’s experiments for more direct comparison. Spheres 
flying at roughly M^ = 14.5 in pure C0 2 for Re^ = 
{ 192, 4564} were simulated with the DSMC code employ- 
ing the same surface description as above. The results are 


included in the drag plots and appear to follow the exper- 
imental trends quite well. Additional simulations will be 
performed to reproduce the Pioneer Venus probe experi- 
ments and results will be reported in the final paper. 

Surface Temperatures and Pyrolysis 

The DSMC calculations above employed a radiative- 
equilibrium energy balance to compute the surface tem- 
perature. This boundary condition neglects any effects of 
heat capacitance and heat conduction in the spacecraft heat 
shield and, therefore, provides an upper bound for the sur- 
face temperature. The calculations also neglected pyrolysis 
from the heat shield which can be significant for a high- 
energy incident flow at very low densities. 

To check the validity and accuracy of these assump- 
tions, the CM A code 27,28 was used to calculate the one- 
dimensional heat transfer into the carbon phenolic heat 
shield at the stagnation point. Temperature -dependent ma- 
terial properties, surface re-radiation, and in-depth pyrol- 
ysis were included in the calculation, but surface ablation 
was neglected. The initial temperature was estimated to be 
150 K just prior to entry from deep space. The aerother- 
mal heat flux was taken from the DSMC results associated 
with A = 0.75 appearing in Table 2. 

Figure 8 presents the transient surface temperature at 
the stagnation point during Galileo entry. As expected, the 
CMA-calculated surface temperature is always below the 
radiative-equilibrium temperature. The temperature differ- 
ence between the two results increases from about 200 K 
at the initial time (Re^ = 0.098) to 1130 K at 55 seconds 
(R^oo = 926). 

Figure 9 compares the stagnation point pyrolysis-gas 
mass flux with the free stream mass flux. The pyrolysis gas 
flux becomes significant between 43 and 48 seconds as the 
material temperature rises from 600 K to 1000 K. At 46.4 
seconds (Re=103.4) the pyrolysis gas flux is more than 20 
times the free stream mass flux. 

Perhaps the drag coefficient does not depend strongly on 
the surface temperature which can itself be significantly be- 
low the radiative equilibrium value during entry. However, 
the surface mass flux due to pyrolysis of the material is sig- 
nificant once the probe drops below roughly 420 km. This 
could lead to a noticeable increase in drag and a decrease in 
heating. Thorough simulation requires that this mass flux 
be coupled directly into the DSMC code. Such modifica- 
tions are underway and the results will be reported in the 
final paper. 

Concluding Remarks 

Drag on the Galileo probe during initial entry into the 
Jovian atmosphere was computed with a DSMC method, 
and the results were consistent with the experimental re- 
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suits. However, drag varies considerably with Reynold’s 
numbers in the range 10 < Re^ < 10, 000 at the high M w 
flight conditions anticipated for the probe. Results corre- 
late better with experiment when plotted against Knudsen 
number rather than Reynolds number under highly rarefied 
conditions. Simple models for gas-gas and gas-surface in- 
teraction permit efficient yet meaningful simulation of the 
flow, although the effects of significant surface pyrolysis 
warrant further investigation. Uncertainties in the thermal 
accommodation coefficient A have minimal impact on final 
results due to insensitivity of drag to A despite great depen- 
dence of heating upon A. Further comparisons to available 
experimental data will be made and the pyrolysis models 
will be coupled directly into the DSMC code before the fi- 
nal paper is presented. 
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Tabic 2: Galileo Probe Entry Simulation Results 


Re 

A 

CD 

q-nose 

T-nose 




(W/m2) 

(K) 

0.10 

0.75 

2.060 

766 

358.0 

1.03 

0.75 

2.048 

6,495 

602.4 

8.69 

0.75 

2.024 

48,728 

998.1 

34.01 

0.75 

1.970 

168,523 

1,365.1 

103.41 

0.75 

1.896 

448,881 

1,737.4 

322.47 

0.75 

1.777 

1,201,262 

2,236.4 

926.15 

0.75 

1.629 

2,980,082 

2,798.3 

0.10 

0.50 

2.061 

474 

318.8 

8.69 

0.50 

2.046 

32,564 

909.9 

103.41 

0.50 

1.934 

309,602 

1,580.3 

926.15 

0.50 

1.628 

2,085,700 

2,545.3 

0.10 

0.90 

2.043 

848 

366.8 

8.69 

0.90 

2.008 

57,941 

1,050.3 

103.41 

0.90 

1.864 

528,677 

1,826.8 

926.15 

0.90 

1.648 

3,559,079 

2,931.1 


Table 1: 

Galileo Probe Entry Simulation Conditions 




t 

(sec) 

Alt. 

(km) 

Re 

M 

Kn 

T 

(K) 

P 

(ksjm2) 

D 

(cells) 

0.00 

735.00 

0.10 

31.47 

415.59 

425.0 

2.024E-11 

16 

19.00 

604.00 

1.03 

34.15 

43.01 

360.5 

1.896E-10 

16 

33.50 

506.00 

8.69 

37.20 

5.55 

303.2 

1.422E-09 

16 

41.00 

453.00 


39.61 

1.51 

267.0 

5.100E-09 

16 


416.00 

103.41 

41.86 

0.53 

238.7 

1.435E-08 

32 

51.25 

382.00 

322.47 

44.58 

0.18 

210.0 

4.097E-08 

32 

55.50 

353.00 

926.15 

46.35 

0.07 

194.0 

1.114E-07 

48 




RolaUonal Degree* o( Freedom 
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Fig. 1 Temperature variation of viscosity of the Jovian atmo- 
spheric gas mixture computed by Biolsi and fit with the VHS 
model. 



Fig. 4 Translational temperature along the stagnation streamline 
ahead of the probe for entry at several values of Re M . 



Fig. 2 Temperature dependence of the rotational degrees of free- 
dom for . 



Fig. 3 Comparisons of probe geometry to simulation configura- 
tions employing differing resolutions. 



Fig. 5 Mass density along the stagnation streamline ahead of the 
probe for entry at several values of Re M . 
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Fig. 6 Drag coefficients from DSMC simulation of Galileo en- 
try compared to results of related experiments of Intrieri; includes 
results from simulations of sphere experiments. Results plotted 
against Re^ . 
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Fig. 8 Radiative-equilibrium and CMA -calculated stagnation 
point material temperature for initial entry period of the Galileo 
probe. 
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Fig. 7 Drag coefficients from DSMC simulation of Galileo en- 
try compared to results of related experiments of Intrieri; includes 
results from simulations of sphere experiments. Results plotted 
against Kn^ . 



Fig. 9 Comparison of frees cream mass flux and CMA -calculated 
pyrolysis gas mass efflux from the probe surface during inital entry 
phase. 





